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Abstract 

We present predictions for /i —> ey and r —» /17 in the CMSSM in which there is 
a natural hierarchy of neutrino masses resulting from the sequential dominance of 
three right-handed neutrinos. We perform a global analysis of this model in the 
(mo, Mi / 2 ) plane, assuming radiative electroweak symmetry breaking, and including 
all observed laboratory data. We confirm that a large (small) r —> fi'y rate results 
from the dominant right-handed neutrino being the heaviest (lightest) one. We show 
that the p —> ey rate may determine the order of the sub-dominant neutrino Yukawa 
couplings in the flavour basis, but may also be sensitive to effects beyond the leading 
log and leading mass insertion approximations. We also show that the n —> ey rate 
is independent of # 13 , but measurement of this angle may determine a ratio of sub¬ 
dominant Yukawa couplings. 
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1 Introduction 


Atmospheric and solar neutrino experiments have provided convincing evidence for 
neutrino masses and mixings. The current minimal best fit interpretation of the data 
involves three light neutrinos with [1] 

|A m 2 32 \ = (1.7 — 3.3)10 _3 eV 2 , sin 2 20 23 = 0.93 - 1.0, (1) 

A m 2 21 = (4 - 10)10 _5 eV 2 , sin 2 2d i2 = 0.71 - 0.89, (2) 

where = m 2 — m 2 are the neutrino mass squared differences and the approximate 
lcr ranges are shown. These masses and mixings represent the first solid evidence of 
new physics beyond the Standard Model. 

The most elegant explanation of small neutrino masses continues to be the see¬ 
saw mechanism [2, 3]. According to the see-saw mechanism, lepton number is broken 
at high energies due to right-handed neutrino Majorana masses, resulting in small 
left-handed neutrino Majorana masses suppressed by the heavy mass scale. A natu¬ 
ral implementation of the see-saw mechanism which can account for a neutrino mass 
hierarchy and a large atmospheric mixing angle is single right-handed neutrino dom¬ 
inance [4], See also [5, 6]. It is also possible to account for a large solar neutrino 
angle in this framework [7]. For example a large solar angle may result from sequen¬ 
tial right-handed neutrino dominance due to the leading sub-dominance of a second 
right-handed neutrino [7], leading to a full neutrino mass hierarchy m\ <C m 2 -C m 3 . 

In charged lepton interactions of the Standard Model the lepton flavour violation 
(LFV) is suppressed by the heavy see-saw mass scale, or equivalently the small left- 
handed Majorana masses, and is practically unobservable. However in low energy 
supersymmetric (SUSY) theories the situation dramatically changes for the better. 
In SUSY the LFV is imprinted on the slcpton mass matrices, resulting in slepton 
masses which are off-diagonal in flavour space. The off-diagonal slepton masses can 
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then induce LFV in loops suppressed only by the SUSY breaking scale, multiplied by 
ratios of the off-diagonal slepton mass to the diagonal slepton mass [9, 10, 11, 12]. 

The experimental prospects for improving the limits or actually measuring LFV 
processes are very promising. The 90% C.L. limits of BR(r —> /xy) < 1.1 x 10 -6 
[13] and BR(/j —» ey) < 1.2 x 10” 11 [14] are particularly stringent in constraining 
SUSY models. In fact, these limits will be lowered in the next 2-3 years as the 
present B factories, inevitably producing tau leptons along with the b quarks, will 
collect 15-20 times more data and as the new [i —» ey experiment at PSI probes 
the branching ratio down to 1CU 14 [15, 16]. The forthcoming experimental results 
are clearly a powerful incentive for theoretical studies of r —> /xy and /x —> ey, and 
this provides an important underlying motivation for the present paper. It should be 
noted however that LFV violation can originate from other effects unrelated to the 
see-saw mechanism, for example colour triplet scalars in GLIT models, or off-diagonal 
slepton masses generated from flavour violating effects in the SUSY breaking sector. 
In any observation of LFV, such effects would need to be disentangled from the effects 
due to the see-saw mechanism in the CMSSM considered here. 

There is a huge literature on LFV in the framework of the CMSSM, where the 
constraint of universal scalar mass mo and trilinear mass Aq implies that the slepton 
mass non-universality originates entirely from RG effects due to right-handed neutri¬ 
nos between the GLIT scale and the lightest right-handed neutrino mass scale [17]. 
There is a much smaller literature concerned with LFV in models based on single 
right-handed neutrino dominance [18, 19]. In our previous paper [18] we made the 
important observation that a large r —> yxy rate results from models in which the 
dominant right-handed neutrino is the heaviest one, and the neutrino Yukawa matrix 
has a lop-sided form with an order unity Yukawa coupling in the 23 position of the 
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matrix 

/ 0 0 0 \ 

Y v « 001. 

Vo 0 l) 

In subsequent papers [19] a bottom-up analysis of various types of right-handed neu¬ 
trino dominance was discussed and r —> p .7 and /i —> ey were both considered in 
terms of coefficients C l3 which parametrise the off-diagonal slepton masses to leading 
log approximation in the mass insertion approach to LFV. Bounds on the coefficients 
Cij were presented as a function of the sneutrino and second gaugino masses me, and 
M 2 , and Cij were related to Yukawa couplings and right-handed neutrino masses for 
each of the different types of right-handed neutrino dominance [19]. 

In the present paper we provide a dedicated analysis of a particular class of right- 
handed neutrino dominance models, namely sequential dominance (SD). The reason 
why we choose to focus on SD is that it provides a particularly elegant application 
of the see-saw mechanism to the LMA MSW solution, and these are the only models 
where a full neutrino mass hierarchy m\ <C m 2 -C m 3 arises naturally. Another 
attractive feature of SD models is that it may become possible to relate the rate for 
/x —> ey directly to the sub-leading Yukawa couplings. We shall discuss the conditions 
under which this may be possible later. There are several examples of models in the 
literature which rely on SD [20, 5, 21, 22], In the realistic models the dominant right- 
handed neutrino is either the heaviest or the lightest. We refer these cases as heavy 
sequential dominance (HSD) or light sequential dominance (LSD), respectively, and 
discuss these cases separately. 

In this paper, then, we perform a global analysis of HSD and LSD models in the 
framework of the CMSSM, presenting our predictions in the (mo, M \/ 2 ) plane, assum¬ 
ing radiative clectroweak symmetry breaking, and including all observed laboratory 
data. We give predictions for branching fractions for r —> /ry and fi —> ey in terms of 
a convenient parametrisation which we introduce for the neutrino Yukawa couplings. 
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Our results confirm that the r —> p-y rate distinguishes HSD from LSD. We further 
show that the yU —> e 7 rate may allow a determination of the sub-dominant neutrino 
Yukawa couplings in the flavour basis, in terms of an expansion parameter A. We 
show that measurement of # 13 may determine a ratio of sub-dominant Yukawa cou¬ 
plings. We discuss quantitative effects which were not present in the leading log and 
leading mass insertion approximations of [19], for example cases where yU —> ey is 
controlled by the 13 slepton mass. 

The remainder of the paper is set out as follows. In section 2 we discuss sequential 
dominance and introduce our parametrisation. We also make some leading log an¬ 
alytic estimates of off-diagonal slepton masses in terms of this parametrisation, and 
discuss the implications for LFV processes. In section 3 we give our full numerical 
results. 


2 Sequential Dominance 

2.1 Brief Review 


In this subsection we give a brief review of SD, including the analytic estimates of 
neutrino masses and mixing angles in terms of Yukawa couplings and right-handed 
neutrino masses, which we assume here to be real. More details can be found in [ 8 ]. 
Note that all the results in this subsection are independent of the mass ordering of 
right-handed neutrino masses, and so apply to both HSD and LSD. 


In the flavour basis where the charged lepton mass matrix is diagonal, the diagonal 
right-handed neutrino mass matrix is written as 


Mrr — 


( X ' 
0 

V 0 


0 

X 

0 



( 3 ) 
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and the neutrino Yukawa matrix without loss of generality as 

/a' a d \ 

Ylr= b> b e . (4) 

\d c f ) 

In order to account for a neutrino mass hierarchy and large neutrino mixing angles 
in a natural way (without fine-tuning) the following SD condition was proposed [7], 

|e 2 |, |/ 2 |, |e/| _ \xy\ _ \x'y'\ 

-y- (5) 

where x,y G {a,b,c} and x',y' G {a',b',c'}. It is further assumed [4] that 


d <C e « /. 


Then it was shown that the neutrino masses are given to leading order in rn 2 jm- 3 by 


n( x 'y' 2 n 

-v 2 ) 

a 2 + ( c 23 b - s 23 c) 2 -2 

x l ’ 2 
e 2 + f 2 2 


where v 2 is a Higgs vacuum expectation value (vev) associated with the (second) 
Higgs doublet that couples to the neutrinos. Note that with SD each neutrino mass 
is generated by a separate right-handed neutrino, and the origin of the neutrino mass 
hierarchy is thus linked to the sequential condition in Eq.5. The neutrino mixing 
angles are given to leading order in m 2 /m 3 by [8], 


tan d 2 


tan O-i 


(c 23 b - s 23 c ) 
a(s 23 b + c 23 c) Y 
(e 2 + f 2 ) X 


Thus, assuming SD, the atmospheric angle is given by a ratio of dominant right- 
handed neutrino couplings associated with Y, the solar angle is given by a ratio of 
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sub-dominant right-handed neutrino couplings associated with X, and #13 is of order 
m 2 /m 3 . The sub-sub-dominant right-handed neutrino X' with couplings a',b',d is 
completely irrelevant for neutrino physics unless m 1 can be measured. In writing the 
equation for # 13 we have assumed for simplicity that 

d <C a(,§23& + C23c)Y/(J e 2 + f 2 X). (9) 


At leading order in a mass insertion approximation [9, 11] the branching fractions 
of LFV processes are given by 

a 3 

BR(/j -»• Iji) an ^2-/(M 2 ,/r,m £; )|m|.J 2 ^ i tan 2 /3 (10) 

where l\ — e, I 2 — /i, I 3 = r, and where the off-diagonal slepton doublet mass squared 
is given in the leading log approximation (LLA) by 

« _(M±A>cv (11) 


where the leading log coefficients are given by 


C 21 = a'l> In -§7 + ahlu ^Y + de In ^ 

n U '1 Mu , 1 1 Mu 1 Mu 

C 32 = bc\n — + bc\n — + ef\n — 


X' X Y 

r< / /! M u , , Mu M v 

C 31 = a c In — + acln — + df In— 


The factors in Eq.10 represent the ratio of the lcptonic partial width to the total 


width, 


sp' — 


T(/j * IjVjVj') 

T(li -> all) 


Clearly C 21 = 1 but £32 is non-zero and must be included for correct comparison with 
the experimental limit on the branching ratio for r —> /ry. This factor is frequently 
forgotten in the theoretical literature. 


We shall focus on C 2 i and 6*32 which correspond to /1 —> ey and r —> /ry. For HSD 
the couplings a', b', d are expected to be smaller than the couplings a, b, c, and we will 
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therefore drop them in both our analytic and numerical analysis. For LSD the primed 
terms are also not relevant for LFV since as discussed later we will set X' = Mu which 
means that this right-handed neutrino is immediately decoupled at the GUT scale, 
and so In XXl = 0. In the case of LSD this is rather a strong assumption and if 
X' < Mu these Yukawa couplings will again become relevant. Dropping the primed 
contributions for both HSD and LSD the relevant coefficients are then given by 


C21 

C32 


ab In 
be In 


M v 

X 

Mu 

X 


+ de In 

+ e .f l n 


M v 

Y 
Mu 

Y 


(14) 


Note that the primed couplings are irrelevant to either the neutrino masses and 
mixing angles or the LFV processes. This applies to both the HSD and the LSD 
cases. 


2.2 A convenient parametrisation 

The above analytic results for the neutrino masses and mixing angles suggest the 
following parametrisation of the Yukawa couplings and right-handed neutrino masses 
in terms of order unity coefficients a l} and A, together with an expansion parameter 
A raised to integer powers, 


a = ar 2 X n f 
b = a 22 X n f 
c = a 32 X n f 
d = a 13 X m f 
e = a 23 f 
X = AX^Y 


( 15 ) 
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where 


A = 


| Am; 


211 


\ | Am; 


2 I 

321 


0.15. 


(16) 


Currently the data prefers A ~ 0.15, with large uncertainty. It is the smallness of this 
value that is the motivation for sequential dominance. 

In terms of the parametrisation in Eq.15, the neutrino mixing angles in Eq.8 
become, 


tan 9 


23 


023 


tan 9 


®12 


12 


9 


022(023 — -S230) 

A 012022(523 + c 23 -r) 


13 


(17) 


A (1 + CI 23 ) 

where for convenience we have also defined the following ratio of sub-dominant Yukawa 
couplings which clearly plays a crucial role in determining # 13 , 


_ £ _ O 32 
b a 2 2 


(18) 


The motivation for the above parametrisation is that it reproduces the large atmo¬ 
spheric and solar mixing angles tan #23 ~ 1, tand 12 ~ 1, with the dimensionless 
parameters being of order unity. 1 

Similarly the neutrino masses in Eq.7 become, in terms of the parametrisation in 
Eq.15, 


m 1 ~ 0 

m 2 ~ (0.12 + 022(02.3 - s 23 r) 2 )^y 

m 3 ~ (1 + 023)^7 (19) 

The ratio of second and third neutrino masses is then m 2 /'m^ ~ A which motivates 

the definition of the expansion parameter in Eq.16. A further motivation for this 

1 Note that with all and A of order unity this is not the most general parametrisation that 

satisfies the conditions in Eqs.5,6 that define this class of models. Also note that as before we have 
assumed for simplicity that d is small in the equation for $ 13 , and the condition for this is now 
to > 2 . 
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parametrisation is that it clearly shows that the neutrino masses and mixing angles 
in Eq.19, 17 are completely independent of n. Although the neutrino masses and 
mixings cannot determine the choice of integer n, which controls the magnitude of 
the Yukawa couplings, the LFV processes are able to do so as we now discuss. 

Using our parametrisations in Eq.15 the coefficients relevant for [i —> ey and 
r —> n'y from Eq.14 are then determined by the products of Yukawa couplings 

ab = ttv 2 (i 22 ^ n 
de = Gq 3 a 23 A m f 2 
be = a 2 22 rX 2 n f 2 

ef = a 23 / 2 (20) 

Then from Eqs.14, 20 we have 

C 2 i = a 12 a 22 A 2jl / 2 ln^ + a 13 a 23 A m / 2 ln^ 

C 32 = a\ 2 rX 2n f 2 In ~ + a 23 f 2 In ^ ( 21 ) 

The coefficients C t -j which determine the strength of LFV processes are clearly gov¬ 
erned by the value of the Yukawa coupling / together with the integers n, m which 
control the magnitudes of the Yukawa couplings. The values of the right-handed 
neutrino masses does not directly influence the coefficients Cij very much, with only 
a mild logarithmic dependence. However the values of the right-handed neutrino 
masses will have an important indirect effect on LFV processes via the values of the 
Yukawa couplings. For example for a fixed m 3 , U 32 oc f 2 oc Y, so the r —> p .7 rate in¬ 
creases as the dominant right-handed neutrino mass Y becomes heavier [18]. We shall 
now discuss the two cases HSD and LSD separately, corresponding to the dominant 
right-handed neutrino of mass Y being the heaviest or the lightest, respectively. 
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2.3 Heavy Sequential Dominance (HSD) 


HSD is defined by the condition that the dominant right-handed neutrino of mass Y 
is the heaviest, with 

Y > X > X' ( 22 ) 

Typically in unified models the heaviest right-handed neutrino is in the 33 position 
and the 33 element of the neutrino Yukawa matrix is equal to the top quark Yukawa 
coupling f — h t (at the high energy scale). The HSD case then leads to the “lop¬ 
sided” form of neutrino Yukawa matrix with an order unity Yukawa coupling in the 
23 position, and subsequently a large expected rate for r —> /ry [18]. For HSD with 
/ = h t , the atmospheric neutrino mass m 3 ~ ^J\Krr^\ then implies that the heaviest 
(dominant) right-handed neutrino has mass 

Y ~ 2 m 2 /m 3 »3x 10 14 GeV. (23) 

where we have used the GUT scale value of the top quark mass rn t which for large 
tan/3 is about 0.7 times its low energy value. For HSD the condition that I « 7 
implies, from Eq.15, that 

n> 1 (24) 

Recall that we also assumed m > 2. The neutrino matrices in Eqs.3,4, in terms of 
the parametrisation in Eq.15, are summarised below for the case of HSD, 



( 

— 

0 

0 \ 


Mgi D = 


0 

HA 2n_1 

0 Y 

(25) 


K 

0 

0 

1/ 



( 

— 

Q 12 A n 

ai 3 X m \ 


v uHSD 

1 LR 


— 

a 22 X n 

023 ht 

(26) 



— 

ra 22 \ n 

1 ) 



where the blanks indicate entries which are irrelevant for both neutrino masses and 
mixing angles and LFV. The neutrino masses and mixing angles are given by Eqs. 17,19, 
independently of n. 
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In this case it is clear from Eq.21 that r —» /iy, which corresponds to C 32 , is large 
and independent of n, being given approximately by 

C 32 = bi -p- (27) 

On the other hand fi —> ey, which corresponds to C 21 , is smaller and more model 
dependent since it depends on the integers n, m , 

C 21 = Oi 2 « 22 A 2 "h 2 in —+ O i3°23^ m / i t bi -p- (28) 

From our numerical results we shall find that the experimental limit on fi —> ey will 
require n > 2 . Since ab ~ A 2 ” while de ~ A m , if 2 n < m then ab will dominate 
over de, while if 2 n > m then de will dominate over ab. In the second case /i —> ey 
is controlled by the 13 neutrino Yukawa coupling d. We emphasise the fact that 
although /i —> ey is model dependent in this framework we understand precisely its 
origin. 

2.4 Light Sequential Dominance (LSD) 

LSD is defined by the condition that the dominant right-handed neutrino of mass Y 
is the lightest, with 

X' > X > Y (29) 

The LSD case then leads a more “symmetrical” form of neutrino Yukawa matrix with 
no large off-diagonal Yukawa couplings, and hence a small expected rate for r —> p.y. 
Such a form of neutrino Yukawa matrix is consistent with an exactly symmetrical or 
anti-symmetrical or mixed symmetry form, or some more general matrix with small 
off-diagonal couplings. LSD is more complicated to discuss since we need to re-order 
the matrices in Eqs.3,4, to put the heaviest right-handed neutrino of mass X' in the 
third column. Then we shall make a similar assumption for LSD that the 33 element 
of the neutrino mass matrix to be equal to the top quark Yukawa coupling c' = h t . 
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For LSD to set the scale of the right-handed neutrino masses, we need to relate the 
dominant right-handed neutrino Yukawa coupling / to c 7 = hf. In particular, we 
extend our parametrisation to include 

/ = 031 AV (30) 

where p > 0, and <231 is an order unity coefficient. For LSD with d = h t , the 
atmospheric neutrino mass m 3 ~ \ Arrif i2 | then implies that the lightest (dominant) 
right-handed neutrino has mass 

Y « a 2 31 X 2p 3 x 10 u GeV. (31) 

with p > 0. The sub-dominant right-handed neutrino mass is heavier than Y and its 
mass is parametrised by 

X = AA 2n_1 Y « Aa 3 1 A 2p+2n_1 3 x 10 14 GeD (32) 


with, 

n < 0 < p. (33) 

In LSD the heaviest right-handed neutrino mass is X' which must be heavy enough 
to satisfy the sequential condition, corresponding to the second inequality in Eq.5, 


,/2 


\xy | 

X X' 


(34) 


From Eqs. 15,32,34 we find 


X' » A " 1 3 x 10 14 Gel/ 


(35) 


independently of p, n. Eq.35 implies that X’ cannot be (much) below the GLIT scale 
Mu. As indicated earlier we shall assume X' = Mu which implies that there can 
be no LFV effects arising from this right-handed neutrino. In order to ensure that 
X <C X ' by comparing Eqs.32 and 35 we also have the condition 


p + n > 0 


(36) 
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The neutrino matrices in Eqs.3,4, in terms of the parametrisation in Eqs. 15,30 
are summarised below for the case of LSD, after re-ordering them to put the heaviest 
right-handed neutrino in the third column, and the lightest (dominant) right-handed 

neutrino in the first column, 

( a§ 1 A 2p 0 0 \ 

0 Aal 1 X 2p+2n ~ 1 0 3 x 10 u GeV (37) 

Vo o -) 

f a l 3 a 31 \ p+m a 12 a 31 X p+n - \ 

«23«3lA P a-22Cl3lX P+n — h t (38) 

V a 3 iA p ra 2 2 a 3 i\ p+n 1 ) 

where the blanks indicate entries which are irrelevant if X' = Mu- However, as 
we mentioned previously, if X' < Mu such Yukawa couplings will again become 
important. Note that n is a non-positive integer, and in realistic models it is nearly 
always negative so that the Yukawa couplings in the first column are smaller than 
the Yukawas in the second column. The advantage of this parametrisation is that it 
provides a unified parametrisation of HSD and LSD, with n being positive for HSD, 
and negative or zero for LSD. In this unified parametrisation the neutrino masses and 
mixing angles are still given by the same formulas as in Eqs. 17,19, idependently of 
n,p,a 3 1 . 

In this case it is clear from Eq.21 that r —» py, which corresponds to C 32 , is now 
much smaller and model dependent since it depends on the integers p, n, 

C 32 = ra 2 22 a 2 3 l X 2 p+ 2 n h 2 t In ^ + a 22 a 2 31 X 2 p h 2 In (39) 

If 7i < 0 then, in the notation of Eq.20, be dominates over ef. Turning to /i —> ey, 
which corresponds to C 2 \ , since n < 0 while m > 2 we have approximately, 

C 21 = a 12 a 22 a 2 31 X 2 p+ 2 n h 2 In ^ (40) 

Ji 

so that, in the notation of Eq.20, ab always dominates over de. By comparing Eqs.39 
to 40, we see that in all cases 

C 32 ~ C 21 ~ A 2 p+ 2 n h 2 In ^ (41) 


M 


LSD 

RR 


Y, 


vLSD 


LR 
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which leads to the LSD prediction that the branching fractions for r —» p'y and 
p —> e 7 should be comparable, and are controlled by the value of p + n. We already 
saw in Eq.36 that X <C X' requires p + n > 0. We shall see in the next section 
that the experimental limit on p — > ey will require p + n > 2. This implies that the 
t —> p 7 rate for LSD is suppressed relative to that for HSD by of order A 8 or more, 
which effectively makes r —> p'y unobservable in this case. In the case that X' < Mu 
and the 23 Yukawa coupling were set equal to the 32 Yukawa coupling we would find 

C 32 ~\P+ n h 2 t (42) 

which would imply that the r —> p'y rate for LSD is suppressed relative to that for 
HSD by of order A 4 or more, which is still at least three orders of magnitude below 
current limits. 

3 Numerical Results 

Our numerical results are based on the top-down global analysis of the CMSSM with 
universal soft scalar mass m 0 , soft gaugino mass M l / 2 , soft trilinear mass A 0 and the 
sign of the Higgs superpotential mass parameter /j which we take to be positive as 
implied by the recent muon anomalous magnetic moment measurement. In a standard 
notation [23] the complete list of the input parameters varied at the unification scale 
Mu ps 2 x 10 16 GeV (which itself is allowed to vary) contains the soft masses m 0 , 
Mi / 2 , Aq, the unified gauge coupling au, the deviation of the QCD gauge coupling a 3 
from au, e 3 = (a 3 — au)/cxu , the deviation of the b, r Yukawa couplings hb, T from the 
(varying) top Yukawa coupling h t , e^ T = (/r biT — h t )/h t , and the five parameters of the 
neutrino sector Y, a i2 , a 22 , «23 arid A defined earlier. The five neutrino parameters 
are only allowed to vary within 20 — 30% (see Eq.15) but typically good fits are 
obtained for values in the smaller range ~ 10% which means that Y and the a’s stay 
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close to 3 x 10 14 GeV and 1, respectively, r = 032/0122 is held fixed in the analysis . 2 /i 
is an input parameter at the low scale which we take as the M z scale. The function 
X 2 = J2( x f MSSM — x^ xp ) 2 /a 2 is evaluated based on the following observables x t : a, 
Ga :i (M z ), M t , rnb(rrib), M r , M z , M w , Ap, BR{b —> sy), a^ EW and the neutrino 
mass differences and mixing angles of Eqs.(l) and (2). The low-scale value of the 
bilinear parameter Bg is directly related to ratio of Higgs vacuum expectation values 
tan/? and the latter is kept fixed in the analysis. Each computed charged fermion 
mass includes the complete 1-loop SUSY threshold correction and correct radiative 
clectroweak symmetry breaking at one loop is required at all points by adopting tight 
c m z an d <Jm w ■ A direct search limit is applied to the mass of each unobserved particle 
including m h 0 > 114 GeV. We shall present our results as contours in the (m 0 , Mi/ 2 ) 
plane, for fixed tan/?, with H 0 and \g\ varying over the plane. More details of the 
global analysis can be found in [23] and [24], 3 

The main constraint of the present study comes from the muon g — 2. Based on 
[25, 26] we require the contribution from new physics to fit 

a N EW = a Jf xp) - af 1 = (34 ± 11) x 10" 10 . (43) 

The main focus of the analysis is neutrino masses and mixing angles which we fit to 
the central atmospheric and LMA MSW values. $13 is predicted in terms of a ratio 
of Yukawa couplings as discussed above and computed from the exact form of the 
MNS matrix obtained at low energy. We emphasise that the neutrino parameter fit 
is within the framework of SD, and in this framework it becomes possible to make 
direct connections between LFV processes and specific Yukawa couplings, subject 

to the conditions discussed in sections 2.3 and 2.4. This would not be possible for 

2 We study cases with r taking on different values as explained below. 

3 We note that in the presented top-down analysis the first and second generation quark and 
charged lepton masses and CKM elements are always fit very well in a separate minimisation with 
small first and second generation yukawa couplings as input at Mjj. The presence of the 3x3 
yukawa matrices Y u , Y d and diagonal Y e renders the whole analysis more complete and consistent 
while, at the same time, the additional flavour structure of Y u and Y d has no effect on the rest of 
the analysis. 
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example if the hierarchy resulted from the tuning of parameters. It is also worth 
noting that we perform a careful RG analysis, adjusting the RG evolution below each 
right-handed neutrino mass threshold, and also perform an exact one-loop calculation 
of lepton flavour violating processes. Thus we have very accurate predictions for 
the t —> /xy and /i —> ey experiments which is increasingly important dne to the 
experimental progress. We shall compare our results to the leading log and mass 
insertion approximations described in the previous sections, and in some cases find 
large discrepancies. 

In Figure 1 we show the results of a global analysis of the HSD model for tan f3 = 
50, n = 2, r = —1, = 0, and the parameters refer to the matrices in Eq.25,26. 

This choice of parameters we call the “baseline” parameter set. Figure 1(a) shows 
the global y 2 fit across the (m 0 ,Mi/ 2 ) plane, with the best fit in the region m 0 ~ 
Mi /2 ~ 500 GeV. The fit deteriorates for large M\/ 2 and mo because of the muon 
anomalous magnetic moment which is shown in the panel (b) and for small M\/ 2 and 
small mo because of b — > sy. In the first case the sparticles are too heavy to generate 
the observed discrepancy (43) while in the latter case they are too light making the 
chargino contribution to the b —> sy effective amplitude too large. 

In Figure 2 we show the predictions of the HSD model with the baseline parameters 
for the branching ratios of r —> /xy and /x —> ey. As expected from Eqs.27 and 28 
the upper panels show a large rate for (a) r —>• /xy and (b) /x —>• ey close to the 
current limits. The prediction for the large r —> /xy rate in panel (a) is quite a 
robust prediction of HSD [18]. As discussed later the predictions are expected to 
scale as tan 2 f3. The /x —> ey rate in panel (b) is sensitive to the choice of n, which 
we have taken to be n = 2. This rate can easily be enhanced by taking smaller n or 
suppressed by taking larger n, as is clear from Eq.28. In the lower panels of Figure 2 
we investigate the accuracy of the leading log approximation (LLA) used to calculate 
the off-diagonal slepton masses in Eq.ll, with the HSD coefficients in Eqs.27,28, as 
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compared to the exact numerical calculation used to generate our numerical results. 
We define the fractional error as follows: 

2 (LLA) 2 

m f — m f 

A = v - hL (44) 

rn~ 

Lij 

Figure 2(c) shows contours of A 32 while panel (d) shows A 2 i, in each case giving 
a measure of the error that would have been incurred had the LLA been used in 
calculating the rates for r —> yuy and /i —> ey. The rates shown in the upper panels 
were, of course, calculated exactly without using the LLA approximation. The LLA 
can induce errors of up to 50% for the lighter SLTSY masses. Note that the errors 
double for the branching ratio. 

The parameter dependence of BR(yU —» ey) in the HSD model is examined in 
Figure 3, where in each case we allow one or two parameters to change from the 
baseline parameter set used in Figures 1 and 2, in order to illustrate a particular 
effect. In Figure 3(a), we suppress the subdominant Yukawa couplings by taking 
n = 4, leaving all the other baseline parameters unchanged. As expected from Eq.28, 
this has the effect of suppressing BR(yU —> ey) by a relative factor of A 8 ~ 2.5 x 10~' 
compared to the previous result, which approximately corresponds to what is observed 
by comparing Figure 3(a) to Figure 2(b). This result demonstrates how sensitive 
BR(/x —> ey) is to the values of the subdominant Yukawa couplings of order X n , with 
the prospect that a measurement of this rate is equivalent to a measurement of these 
Yukawa couplings parametrised by n. This is an exciting prospect offered by the SD 
class of models, however this result is subject to the discussion below Eq.28, which 
we shall return after the following paragraph. 

We now discuss the numerical importance of effects beyond the leading mass 
insertion approximation. In Figure 3(b) we set a 2 2 = 0 with the other parameters as 
in the baseline set, in particular 013 = 0. Having a 2 2 = a 13 = 0 should completely kill 
/X —» ey according to Eq.28, however clearly it does not, since the rate in Figure 3(b) 
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is in fact larger than the rate in Figure 3(a) by three orders of magnitude! Clearly 
the mass insertion approximation, on which our analytic expections of the previous 
section, and all the results in [19] for example, are based, is no longer reliable in the 
limit of a 2 2 = Ui 3 = 0. We note that having a 22 = Ui 3 = 0 is completely reasonable 
from the SD point of view of neutrino masses and mixing angles, and would correspond 
to some specific type of texture of the “Fritzsch” type [27] for example. Although in 
this case C 2 i = 0, the values of 6*32 and C 31 are non-zero, leading to off-diagonal 32 
and 31 slcpton masses, where according to Eq.12 the 31 slepton mass must be arising 
from the ac product of Yukawa couplings which would be suppressed by increasing n. 
A 21 flavour violation may be generated by a double mass insertion involving both 
32 and 31, and this results in the observed rate for 11 —> ey in panel (b) of Figure 3, 
although we calculate the rate exactly without relying on any mass insertion methods 
at all as in [12, 18], for example. We call this the 13 slepton effect. 

We now consider the effects of the 13 Yukawa coupling which contribute to 
BR(/x —> ey) according to the second term of Eq.28. Such effects are illustrated 
in Figure 3(c) where we keep n — 4 as in panel (a), but now take 013 to be of order 
one, with the 13 Yukawa coupling controlled by the value of m, which we take to be 
m = 3. Switching on the 13 Yukawa coupling in Figure 3(c) results in a dramatic 
increase in BR,(/i — > ey) compared to Figure 3(a). We call this the 13 Yukawa effect. 

Figure 3(d) shows BR(yU —> ey) in the HSD model where tan/3 = 10 and the 
remaining parameters take the baseline values. Comparing to the baseline contours 
in Figure 1(b), it is clear that the contours do not simply scale as tan 2 /?. The reason 
is that \n | is fixed by electroweak symmetry breaking, and also depends on tan/3, 
and this leads to the difference in the two sets of contours. If // were held fixed 
then we would find the expected tan 2 /? scaling of the rates. The figures illustrate 
the importance of the implicit \n\ dependence for differing values of tan/3, and warn 
against naively assuming a tan 2 /? scaling of the results in the (mo,M!/ 2 ) plane. 
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In Figure 4 we show the prediction of sin 2 2ffi 3 as a function of the ratio of sub¬ 
dominant Yukawa couplings r = 032/022 for the HSD class of models, taking all the 
other parameters equal to their baseline values, in particular 033 = 0. The qualita¬ 
tive variation of this angle with r follows from Eq.17 which shows that although the 
natural value of the angle 6 1 3 is set by the parameter A, as r —> —1, 9 13 —> 0, while 
maintaining tan #12 ~ 1. Numerically we perform a global fit for three specific values 
of r and for each value we find a band of values of sin 2 2 #i 3 corresponding to the vari¬ 
ation in the values of the order one coefficients. For r — — 1 the value of sin 2 20i 3 is 
not exactly equal to zero, due to subleading effects beyond the level of approximation 
of Eq.17, but it does get very small. In general it is fair to say that over most of 
the range of r the value of sin 2 20 13 is below the current CHOOZ limit of about 0.14 
but for r > 0 it should be within the combined ICARUS and OPERA expected 2009 
limit of 0.03. However it is clear that for r < 0 a measurement of sin 2 20i 3 would 
only be possible at JHF or a Neutrino Factory. Note that Figure 4 assumes Gq 3 = 0, 
and the result is more complicated if this is not the case. We emphasise that for 
«i 3 = 0 the rate for /i —> ey is generally independent of the prediction for 6 * 13 , which 
is clear analytically by comparing Eq.17 to Eq.28. The point is simply that C' 2 i is 
independent of a 3 2 and hence r. 4 We emphasise that #13 can be rather large while 
/i —> ey can be arbitrarily suppressed simply by increasing the value of n. 

In Figure 5 we show predictions of the LSD model for the branching ratios of 
r —> //y and ft —> ey. The parameters used are tan fj — 50, p + n — 2, n = —1, 
r = — 1 , 013 = 0, where the parameters refer to Eqs.37,38. As before, the remaining 
parameters A, aij are allowed to vary over the range 0.7—1.3 (and again typically good 
fits are obtained for values in the smaller range 0.9 — 1.1). The LSD prediction for 

BR(/x —> ey) in Figure 5(b) is comparable to that for the HSD model in Figure 2(b). 

4 In the finite 13 Yukawa case 813 is given by a more complicated formula [8]. However, the 
experimental limit on — > ey) is a strong constraint in this case. 
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This is expected since the coefficient C 21 for the LSD model in Eq.40 with p + n — 2 
is of similar magnitude to the coefficient C '21 for the HSD model in Eq.28 with n = 2. 
The new feature in the LSD model is that the magnitude of BR(r —> p'y) in Figure 
5(a) is now much smaller, and in fact is similar to the magnitude of BR(yU —> eq) 
in Figure 5(b), as predicted by Eq.41. However, if we had allowed X’ < Mu then a 
larger, but still suppressed, prediction for r —» p'y following from Eq.42 would have 
been found. In any case the practical conclusion is the same: in the forseeable future 
r —» pq is experimentally unobservable for the LSD case. However p — » eq may 
be observable and is determined by the choice of p + n, which we have taken to be 
p + n = 2. Figures 5(c),(d) show the accuracy of the LLA for the LSD class of models. 
As before the accuracy improves for larger SUSY masses but the error doubles in the 
evaluation of the branching ratio. 

Note that in the case of the LSD class of models the branching ratio for p —> eq can 
easily be enhanced by taking smaller p + n or suppressed by taking larger p + n, as is 
clear from Eq.40. Thus, a measurement of this branching ratio gives a measurement of 
the subleading Yukawa couplings of the second column, as in the HSD case. However, 
unlike the HSD case, the LSD case is not sensitive to the 13 slepton effect and the 
13 Yukawa effect in Figures 3(b),(c). The effects beyond the leading mass insertion, 
due to a double slepton mass insertion, and effects coming from the Yukawa coupling 
parametrised by 031 , are both negligible for the LSD case. The reason that these 
effects were important for the HSD case was due to the large, order unity, Yukawa 
coupling in the 23 position. In the LSD case there are no large Yukawa couplings 
apart from that in the 33 position, and hence such effects are not so important in 
this case. This means that in the LSD case, the measurement of BR(yU —> eq) is more 
reliably related to a measurement of the Yukawa couplings parametrised by p + n. 
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4 Summary 


This paper has focused on models which satisfy the sequential dominance condition. 
We have considered separately two broad classes of models, HSD and LSD, in which 
the dominant right-handed neutrino is the heaviest or lightest one. Although specific 
examples of models exist in the literature which satisfy either the HSD [21] or LSD [22] 
conditions, we emphasise that the sequential dominance condition is simply equivalent 
to assuming that a neutrino mass hierarchy is generated via the see-saw mechanism 
without fine-tuning in the leptonic sector, in the presence of large atmospheric and 
solar mixing angles. Therefore our results in fact apply to a very large class of models, 
although some of our specific results apply only in certain specific regions of parameter 
space. 

We have presented approximate analytic results for neutrino masses and mixing 
angles, and for off-diagonal slepton masses in the leading log approximation in each 
case, using a new parametrisation of the matrices at the high energy scale. The 
parametrisation is very useful in relating the neutrino and LFV observables and in 
interpreting the data in terms of model parameters. 

The numerical results, extracted from the best fits of the global top-down anal¬ 
ysis, demonstrate the validity (limited, at times) of the leading log and leading 
mass insertion approximations used in the literature. We find that the magnitude 
of BR(r —> frf) provides the main discriminator between the different classes of se¬ 
quential dominance, with a large rate not far below current limits predicted by HSD, 
and a much smaller rate predicted by LSD. The observation of BR(/x —> ey), on the 
other hand, may determine the order of the sub-dominant neutrino Yukawa couplings 
in the flavour basis in both the HSD and LSD cases. However, due to the large 23 
Yukawa coupling, the results for HSD may be quite different from the leading log 
and leading mass insertion approximation predictions, making the interpretation of 
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BR^yx —> ey) in terms of underlying Yukawa couplings more difficult, but still pos¬ 
sible in principle. We have also shown that BR^/x —> ey) is independent of $ 13 , but 
measurement of this angle may determine a ratio of sub-dominant Yukawa couplings. 

Ultimately one wants to use the observed neutrino data and constraints (or, if we 
are lucky, positive results) from the measurements of the BR(r —»• /xy) and BR(/x —> 
ey) and other LFV processes to construct models and discriminate among them. To 
this end our study will serve as a framework for such top-down probes, and should 
help to establish: (i) if a sequential dominance mechanism is at work in the neutrino 
sector; (ii) if so then what type of dominance (HSD or LSD) is relevant; (iii) and then 
(subject to the exceptions discussed in section 3) the magnitude of the sub-leading 
neutrino Yukawa couplings in the flavour basis. 

We emphasise how, within the framework of sequential dominance, BR(r —> /xy) 
and BR,(/i —> ey) have a direct interpretation in terms of the underlying Yukawa 
couplings. It is worthwhile reiterating, however, that if there are additional sources 
of flavour violation, as generically expected in string theory [28], then these effects 
will need to be disentangled from the effects due to the RG running of the see-saw 
mechanism in the CMSSM discused here. Also for simplicity we have neglected the 
effects of phases in a first global analysis of the CMSSM with a natural neutrino mass 
hierarchy, although our results demonstrate that the sign of Yukawa couplings can 
play an important part in determining the mixing angles. 
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m 0 [GeV] m 0 [GeV] 


Figure 1: Results for the global analysis of HSD. We use the baseline parameters 
tan/3 = 50, n = 2, r = —1, = 0. Panel (a) shows the \ 2 of the fit which is 

minimised for m 0 ~ M\i 2 ~ 500 GeV. Panel (b) shows the SUSY contribution to the 
anomalous magnetic moment of the muon. 
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Figure 2: The upper panels show the predictions for the branching fraction for (a) 
r —> fi'y and (b) /i —> e 7 for HSD using the same baseline parameters as in Figure 1. 
The lower panels (c), (d) show the fractional error Ay, defined in the text, that would 
be made in calculating the off-diagonal slepton masses if the leading log approximation 
had been used instead of the exact calculation. 
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m 0 [GeV] m 0 [GeV] 

Figure 3: Predictions for BR(/i —> e 7 ) in HSD for alternative choices of parameters 
in this class of models. In panel (a), we choose tan/3 = 50, n = 4, r = —1, a i3 = 0. 
In panel (b) we choose tan/3 = 50, n = 2, r = 0, ai 3 = 0, a 2 2 = 0. In panel (c) we 
choose tan/3 = 50, n = 4, r = —1, m = 3. In panel (d) we choose tan/3 = 10, n = 2, 
r = - 1 , a 13 = 0. 
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Figure 4: The prediction of sin 2 20 13 as a function of the ratio of subdominant Yukawa 
couplings r = 032/022 for the HSD class of models. The remaining parameters are as 
in Figure 1. 


29 






m 0 [GeV] m 0 [GeV] 

Figure 5: The upper panels show the predictions for the branching fraction for (a) 
r —> /ny and (b) n —► e'y for the LSD class of models. The parameters used are 
tan/? = 50, p + n — 2, n = —1, r = —1, ai 3 = 0. The lower panels (c), (d) show 
the fractional error defined in the text, that would be made in calculating the 
off-diagonal slepton masses if the leading log approximation had been used instead of 
the exact calculation. 
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